Unclassified 

SFCuHlTY 'CLASSIFICATION  OP  FHlS  PAGE 


1.1  K£  POK  T  SECURITY  CLASSIFICATION*^ 

u-cla5S,1.f.t!d _ sgfe  El 

2 a.  SECURITY  ClASSlFlCATlOftJikUlilORlTr 


_  AD-A221  068 

tD'9Q0MENTAT10N  PAGE  I  o 


! 


m  Approved 
OMB  No.  0704-0138 


lb  RESTRICTIVE  MARKINGS 


EC1 E 


2b.  DECLASSIFICATION /OOV 


jING  SCHEDULE 


3  .  DISTRIBUTION /AVAILABILITY  OF  REPORT 

Approved  for  public  release;  Distribution 
unlimited 

5.  MONITORING  ORGANIZATION  REPORT  NUM8ER(S) 


4.  PERFORMING  ORGANIZATI  JWtf  PORT  NU/|aSqS) 5  MONITORING  ORGANIZATION  REPORT  NUI 

CL-TR-Q0-0097  U  C6 

ba.  NAME  OF  PERFORMING  ORGANIZATION  |  bb  OFFICE  SYMUOl  7a.  NAME  OF  MONITORING  ORGANIZATION 

(If  applicable) 

Geophysics  Laboratory  OPB 


I  Geophysics  Laboratory 

[  be.  ADDRESS  (Cify,  State,  and  ZIP  Code) 

Han scorn  AFB 

i  Massachusetts  01731-5000 


7b.  ADDRESS  (C/ty,  State,  and  Z/PCodeJ 


R.i  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 


8b  OFFICE  SYMBOL  9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(If  applicable) 


Be.  ADDRESS  (City,  Stale,  and  ZIP  Code) 


10.  SOURCE  OF  FUNDING  NUMBERS 

PROGRAM  I  PROJECT  1  TASK 

ELEMENT  NO  NO.  NO 

63220C  S321  13 


WORK  UNIT 
ACCESSION  NO. 

01 


ty0 FJ S HAR&n,J  THE  STRATEGIC  HIGH  ALTITUDE  RADIANCE  CODE,  AND  DESCRIPTION 
j  OF  ITS  NEW  AURORAL  MODULE  _ 

I  feilgf 'H5S^tMl)TjRaS6gy!M/,«laliggoggIf*1inl8l?ience  Bernsteln*’  John  Grunlnf.er*, 

j  Ijj,  TYPE  Of  REPORT  13b.  TIME  COVEREO  14.  DATE  OF  REPORT  (Year,  Mon f/l,  Djy)  15.  PAGE  COUNT 

Reprint _  from  to  1990  April  20 _  13 _ 

i6.  supplementary  notation  *Spectral  Sciences,  Inc.  It  Yap  Analytics,  Inc. 

Reprinted  from  Presentation  at  an  IRIS  Meeting  in  Monterey,  CA,  January  1990 

T7  COSATI  CODES  . .  18.  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

- - group - sua ,.(iRQ.u_p -  SHARC;  Infrared  Radiation;  Nonequilibrium  Processes 


17. 

COSATI  CODES  ] 

FIELD 

GROUP 

SUB-GROUP 

— ,.b  The\£tr*cegic\High-Altitude\Radiar.ce\Code  (SHARC)  is  a  new 
computer  code  that  calculates  atmospheric  radiation  in  the 
meaoahpere  and  thermosphere.  The  initial  version,  SHARC-1,  ia 
available  for  distribution.  This  talk  discusses  the  capabilities 
of  this  code  and  describes  the  new  auroral  model  which  will  be 
incorporated  in  the  next  version.  SHARC  calculates  radiance  and 
transmittance  for  paths  from  60  to  300  km  altitude  in  the  2-40  Mm 
spectral  region.  It  models  radiation  due  to  NLTE  (Non-Local 
Thermodynamic  Equilibrium)  molecular  emissions  which  are  the 

dominant  sources  at  these  altitudes.'  Yx’'  .  . 

i  .,irU-a  ,  AY-'p-cf-  M- 


20.  DISTRIBUTION /AVAILABILITY  OF  ABSTRACT 
□  UNCLASSIFIED/UNLIMITED  13  SAME  AS  RPT. 
22a  NAME  Of  RESPONSIBLE  INDIVIDUAL 

Ramesh  Sharma 
DO  Form  1473,  JUN  86 


~.c ,  a 


21.  ABSTRACT  SECURITY  CLASSIFICATION 

□  otic  users _ Unclassified  _ 

22b  TELEPHONE  (Include  Area  Code)  22c.  OFFICE  SYMBOL 

_ (617)  377-4718  OPR 

Previous  editions  are  obsolete.  SECURITY  CLASSIFICATION  OF  THIS  PAGE 

Unclassified 


90  05  02-  019 


DISCLAIME1 NOTIC 


QUALITY  AVAILABLE.  THE  COPY 
FURNISHED  TO  DTIC  CONTAINED 
A  SIGNIFICANT  NUMBER  OF 
PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY. 


GL-TR-90-0097 


CURRENT  STATUS  OF  SHARC, 

THE  STRATEGIC  HIGH  ALTITUDE  RADIANCE  CODE, 

AND 

DESCRIPTION  OF  ITS  NEW  AURORAL  MODULE 
January  1990, 

Raraesh  Sharma,  Anthony  Ratkowski 
Geophysics  Labor atory/OPB , 

James  Duff,  Lawrence  Bernstein,  John  Gruninger, 
Robert  Sundberg  and  David  Robertson, 
Spectral  Sciences,  Inc., 
and 

Rebecca  Healey, 

Yap  Analytics,  Inc. 


•K'C  ■. 


4r 


;  nhs  cka*i 
i  OTIC  TAB  □ 

Ju'JtitiC.iUon  __ _ 


By 

j 

Distribution  | ; 

,  .  _ 

Availability  Codes 

•  — . — 

* 

Avail  and  1  o:  ! 

Dlst 

Special 

‘ 

M-l 

ABSTRACT 

The  Strategic  High-Altitude  Radiance  Code  (SHARC)  is  a  new 
computer  — code  that"  calculates  atmospheric  radiation  in  the 
mesoshpere  and  thermosphere.  The  initial  version,  SHARC- 1,  is 
available  for  distribution.  This  talk  discusses  the  capabilities 
of  this  code  and  describes  the  new  auroral  model  which  will  be 
incorporated  in  the  next  version.  SHARC  calculates  radiance  and 
transmittance  for  paths  from  60  to  300  km  altitude  in  the  2-40  Mm 
spectral  region.  It  models  radiation  due  to  NLTE  (Non-Local 
Thermodynamic  Equilibrium)  molecular  emissions  which  are  the 
dominant  sources  at  these  altitudes. 


1 .  INTRODUCTION 


SHARC,  the  Strategic  High-Altitude  Atmospheric  Radiance  Code, 
calculates  background  IR  radiation  at  high  altitudes. ^ it  includes 

IT  tC.  Sharma ,  A.  Ratkowski,  P.  Acharya,  L.  Bernstein,  J.  Duff,  J. 

Gruninger,  D.  Robertson,  and  R.  Sundberg,  "SHARC  -  An  Atmospheric  Model 
for  High  Altitudes  (60-300  km)".  Presented  at  the  1989  Meeting  of  the 
IRIS  Targets,  Backgrounds  and  Discrimination  (February  1989). 

2.  R.  Sharma,  A.  Ratkowski,  R.  Sundberg,  J.  Duff,  L.  Bernstein,  P. 

Acharya,  J.  Gruninger,  and  D.  Robertson,  "Description  of  SHARC,  The 
Strategic  High-Altitude  Radiation  Code,"  Rpt.  No.  GL-TR-89-0229 , 
Geophysics  Laboratory,  Hanscom  AFB,  MA  01731  (August  1989). 


observer  1 ine-of — s ight  (LOS)  paths  above  60  kra  and  calculates  radiation 
from  the  five  most  important  IR  molecules,  NO,  CO,  H20,  O3  &  C02 .  This 
paper  reports  on  work  done  on  the  model  over  the  last  twelve  months.  The 
initial  version,  SHARC-1,  has  been  made  available  to  the  DOD 
community.1-2)  The  major  upgrade  has  been  development  of  an  auroral  module, 
which  is  based  on  the  GL  AARC  (Auroral  Atmospheric  Radiation  Code).^)  It 
predicts  enhancement  due  to  electron  dosing  during  auroral  conditions.  In 
the  SIIARC  formulation  radiation  from  the  quiescent  atmosphere  around  the 
aurora  is  also  included. 

At  these  altitudes,  most  atmospheric  IR  radiation  comes  from  NLTE 
(Non-Local  Thermodynamic  Equilibrium)  emissions  of  molecules  in  excited 

(i 

vibrational  states.  '  The  NLTE  conditions  arise  because  higher  vibrational 

1  i 

states  that  result  from  chemical  reactions  in  the  upper  atmosphere  do  not 
undergo  enough  collisions  to  restore  thermal  equilibrium  before  they 
radiate,  e.g.,  their  radiative  lifetimes  are  less  than  or  equal  to  their 
colliuional  lifetimes.  Thus,1  the  populations  of  each  vibrational  state 
must  be  calculated  separately.  Figure  1  show  two  calculations  for  the 
spectral  radiance  from  the  15  pm  C02(v2)  band  for  a  space-based  observer 
looking  through  the  earth  limb  with  a  tangent  altitude  of  150  km.  The  NLTE 
radiation  is  much  less  than  that  from  assuming  (incorrectly)  that  the  band 
is  in  LTE .  Referring  to  temperatures  instead  of  vibrational-state 
populations,  the  C02  band  is  vibrationally  much  colder  than  its  kinetic 
temperature.  The  rotational  states  within  each  vibrational  band  are 
assumed  to  be  in  LTE;^'^  this  has  interesting  ramifications  which  will  he 
discussed  below. 

j"!  Winick,  R.  Picard,  R.  Joseph,  R.  Sharma,  and  P.  Wintersteiner , 
"AARC:  The  Auroral  Atmospheric  Radiance  Code,"  Rpt.  No. 

AFGL-TR-87-0334 ,  Air  Force  Geophysics  Laboratory/OPE,  JIanscom  AFB,  MA 
01731  (November  1987)  ADA202432. 

4.  Handbook  of  Geophysics  and  the  Space  Environment,  A.  S.  Jursa  Ed.,  Air 

Force  Geophysics  Laboratory,  Hanscom  AFB,  MA  01731  (1985);  copies 

obtained  from  NTIS,  5285  Port  Royal  Road,  Springfield,  VA  22161,  ADA 
167000. 

5.  R.  Sharma,  "Infrared  Airglow,"  Progress  in  Atmospheric  Physics,  Eds. 
R.  Rodrigo  et  al.,  Kluwcr  Academic  Publishers  (1988). 
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150  km  LIMB  RRDIRNCE 


FREQUENCY  (on'1) 


Figure  1.  Calculation  of  the  Spectral  Radiance  from  the  15  Urn 
C'°2^v2^  Band  Under  Assumed  LTE  (Higher)  and  NLTE 
(Conditions).  The  NLTE  Calculation  is  the  Correct  One, 


2.  SHARC-1  STATUS 

SHARC-1  was  described  at  the  1989  IRIS  TBD  Meeting.  This  is  the 

version  which  is  now  available  to  authorized  users.  Some  of  its  basic 
features  are: 


60-300  km  altitude  regime, 

2-40  pm  with  a  spectral  resolution  of  0.1  cm~*, 

Adjustable  chemical  kinetic  mechanisms  and  rates, 

User-selected  molecular  bands,  and 
Arbitrary  viewing  geometries. 

Arbitrary  LOS  paths  are  allowed,  so  long  as  they  are  above  60  km.  The  next 
version,  scheduled  for  release  at  the  end  of  the  year,  will  extend  down  to 
50  km.  The  data  bases  for  the  atmospheric  profiles  and  the  chemical 
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SHARC 


kineci.cs  are  based  on  chose  of  the  predecessor  code  HAIRM.^*^^ 
consists  of  six  basic  modules: 

INPUT: 

Menu  driven 

Inceraccive  or  batch  mode 


GEOMETRY  (GEOM) : 

Calculates  observer  LOS  geometries 

CHEMICAL  KINETICS  (CHEMKIN) : 

Solves  steady  state  rate  equations 

I 

RADIATION  TRANSFER  (NEMESIS): 

Radiative  excitation  via  layer-layer  coupling 
Radiative  trapping 

SPECTRAL  RADIATION  TRANSPORT  (SPCRAD): 

Single  line  NLTE  equivalent-width  formalism 
-  HITRAN  line  file  augmented  with  NLTE  03  lines 


OUTPUT: 

Excited  vibrational  state  populations  &  temperatures 
LOS  Spectral  and  in-band  radiance 


An  illustrative  SHARC-1  calculation  for  a  limb-viewing  path  with  a 
tangent  altitude  of  80  kra  is  shown  in  Fig.  2.  The  strongest  radiators  are 
the  molecular  bands  of  4,3pm  (^(v^),  15pm  003^2)  and  9.6pm  03(^3).  Other 

6.  T.  Degges  and  A.  D'Agati,  "A  User's  Guide  to  the  AFGL/  Visidyne  High 

Altitude  Infrared  Radiance  Model  Computer  Program,"  Rpt.  No. 

AFGL-TR-85-00 1 5 ,  Air  Force  Geophysics  Laboratory,  Hanscom  AFB,  MA  01731 
(1984)  ADA161432 . 

7.  R.  Sundberg,  D.  Robertson,  R.  Sharma,  and  A.  Ratkowski,  "HAIRM-87  A 
High  Altitude  Infrared  Radiance  Model,"  Rpt.  No.  AFGL-TR-88-0014,  Air 
Force  Geophysics  Laboratory,  Hanscom  AFD,  MA  01731  (1988). 
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LIMB  VIEW  RT  80  km 


WAVENUMBER  (cm-1) 


Figure  2.  Atmospheric  Spectral  Radiance  for  a  Limb  Viewing  Path  with 
a  Tangent  Altitude  of  80  km  Under  Daytime  Conditions. 

molecular  bands  are  indicated  in  the  figure.  The  C02(v3)  band  shows  a 
strong  diurnal  effect  with  the  daytime  emissions  being  an  order  of 
magnitude  stronger.  This  arises  from  solar  excitation  in  the  v^+v^ 
combination  band  at  2.7  urn  and  subsequent  relaxation  via  energy  transfer  to 
the  band.  The  03^3)  band  also  shows  diurnal  structure,  but  it  is 
stronger  at  night.  In  this  case  solar  radiation  breaks  up  the  weakly  bound 
O3  to  form  02  and  0.  It  reforms  during  the  night.  The  radiation  in,  the 
16-30  pm  spectral  region  is  from  rotational  bands  of  ground-state  H20. 


8.  For  example,  see  R.  Sharma  and  P.  Winter  Steiner ,  "C02  Component  of 

Daytime  Earth  Limb  Emission  at  2.7  Micrometers,"  J.  Geophys.  Res., 
90,  9789-9803  (1985),  and  references  therein.  -  —  ~ 


1976  STANDARD  ATMOSPHERE 


Figure  3.  Atmospheric  spectral  Radiance  for  Limb  Viewing  Paths  with 
Minimum  Altitudes  of  120,  160,  and  200  km. 

Three  calculations  for  limb  viewing  with  minimum  LOS  path  altitudes  of 

f 

120,  160  and  200  km  are  shown  in  Fig.  3.  This  figure  illustrates  some 
problems  encountered  in  developing  quantitative  models  for  these  high 
altitudes.  The  two  strongest  radiators  in  this  spectral  region  are  the 
C02(v2>  band  and  rotational  H2O.  The  H2O  emissions  are  much  weaker  than 
those  from  CO2  except  at  the  highest  altitude,  where  they  are  almost  as 
strong.  Determining  whether  this  is  a  model  artifact  requires  further 
study.  As  illustrated  in  Fig.  1  the  CO2  band  is  strongly  NLTE,  while 
rotational  H2O  is  modeled  as  LTE.  It  is  generally  estimated  that  there  are 
sufficient  collisions  up  to  about  150  km  altitude  to  maintain  rotational 
LTE  within  each  molecular  vibrational  state. There  are  neither  data 
nor  models  currently  available  to  describe  at  what  altitudes  and  in  what 
manner  the  rotational  transitions  depart  from  LTE.  The  rapidly  increasing 
ambient  temperature  compensates  for  the  decreasing  H2O  number  density  and 
leads  to  only  a  slow  decrease  in  radiation  with  altitude.  The  NLTE  CC^Cv^) 
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radiation  decreases  roughly  with  the  product  of  the  0  and  its  number 
densities.  At  200  km  the  CO2  radiation  has  decreased  to  the  same  intensity 
as  the  II2 0  rotational  bands,  if  the  latter  are  reasonably  described  by  the 
LTE  model.  This  figure  illustrates  some  of  the  uncertainties  and  unknowns 
involved  in  predictive  models  at  these  high  altitudes;  measurements  are 
very  difficult,  hence  sparse  and  frequently  questionable. 

3.  THE  AURORA  MODEL 

The  auroral  module  calculates  the  enhanced  IR  radiation  from  NO,  CO2 
(4.3  Mm)  and  NO4,,  as  well  as  contributions  from  unexcited  ambient 
molecules.  The  module  has  the  same  user-friendly  structure  as  SHARC-1  in 
that  it  has  an  interactive  input  module  and  all  reactions  and  reaction 
rates  are  accessible  to  the  user.  The  calculation  begins  with 
electron-deposition  models  for  Class  II,  III  and  III+  auroras  of  finite 
size.  The  scientific  approach  used  in  AURORA  is  substantially  based  on 
AARC. 

The  AURORA  module  is  currently  just  being  completed  and  works  in 
conjunction  with  the  other  code  modules.  The  AURORA  module: 

•  calculates  the  deposition  rate  of  the  primary  electrons, 

•  solves  the  time  and  energy  dependent  rate  equations  for  the 
secondary  electron  distribution, 

•  solves  the  chemical  rate  equations  for  vibrational  excitation 
of  C02  and  the  formation  of  vibrat iona 1 ly  excited  NO  and  N0+, 
and  then 

•  passes  the  vibrational  state  populations  on  to  SPCRAD  for 
calculating  the  spectral  radiance  for  the  LOS  passing  through 
the  aurora. 

The  LOS  information,  which  includes  state  populations  and  column 
densities  for  both  auroral  and  ambient  regions  along  the  LOS,  required  a 
reformulated  geometry  module.  Species  populations  for  the  undisturbed 
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background  are  calculated  in  the  usual  manner.  Ambient  and  auroral 
populations  are  then  stored  for  each  layer  along  the  LOS  for  use  by  the 
radiative  transport  module  (SPCRAD)  to  calculate  the  path  radiance  and 
transmittance. 

The  time  integration  of  electron  dosing  uses  the  Gear  method  for  stiff 
differential  equations,  which  is  part  of  the  Sandia  CI1EMKIN  code^9^ 
incorporated  in  SHARC.  The  energy  deposition  model  for  primary  electrons 
is  based  on  work  by  Griin,  Rees,  and  Strickland. (10-12)  an  exampie  Qf  the 
primary  electron  energy  deposition  model,  the  total  ion  pair  production 
rate  for  a  class  III+  aurora  is  shown  in  Fig.  4.  The  actual  routines  were 
taken  from  AARC  and  modified  to  be  consistent  with  the  SHARC  structure.  In 
addition  to  the  time-dependent  auroral  electron  spectra,  AURORA  also 
calculates  the  time  dependent  production  rates  for  the  three  auroral 
species.  The  chemical  reactions  and  energy  transfer  processes  used  in 
these  calculations  are  described  in  the  AARC  manual. 

From  the  ion  nair  production  rate  as  a  function  of  altitude, 
ionization  efficiences,  secondary  electron  distributions  and  electronic  and 
vibrational  excitation  rates  for  the  major  atmospheric  constituents  are 
calculated.  An  illustrative  secondary  electron  distribution  is  shown  in 
Fig.  5.  An  interesting  feature  is  the  large  dip  around  3  eV,  which  is  due 
to  secondary  electrons  transferring  their  energy  to  form  vibrationally 

9.  R.  Kee,  J.  Miller,  and  T.  Jefferson,  "CHEMKIN:  Problem- 
Independent,  Transportable,  Fortran  Chemical  Kinetics  Code  Package," 
Sandia  Rpt.  No.  SAND80-8003,  Sandia  National  Laboratory,  Livermore, 

CA  94550  (March  1980) . 

10.  A.  Griin,  "Luminescenz-photometrische  Messungen  der 
Energieabsorption  im  Strahlungsfeld  von  Elektronenquel len 
Eindimensionaler  Fall  im  Luft,"  Z.  Natiirforsch. ,  112a,  89  (1957) 

11.  M.  Rees,  "Auroral  Ionization  and  Excitation  by  Incident  Energetic 
Electrons,"  Planet.  Space  Sci.,  Ill,  1209  (1964). 

12.  D.  Strickland,  D.  Book,  T.  Coffey,  and  J.  Fedder, 

"Transport  Equation  Techniques  for  the  Deposition  of  Auroral 
Electrons,"  J.  Geophys.  Res.,  181,  2755  (1976). 
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SHARC  (CL)  NO  EFFICIENCIES  <!3C  III) 


o  2  «.  6  e  ia  i?  i« 

Efficiency  (phol ons/ i on -po i r) 


figure  6  Comparison  of  Time  Oepehdent  SHARC  «m*  Steady  State  AAKC 
E f f t < t ene lea  for  NO. 

exerted  Mj,  The  f««tu>  ft  and  altitude  dependence  of  secondary  electron 
dpitrtliutiuni  have  been  discussed  in  detail  e  !  sewher  e  .  **  ^  ) 

Figures  ft  and  7  compare  SHARC  and  AARC  steady  state  production 
cfficiences  for  HO  and  NO*  in  a  Class  III  aurora.  The  calculated  time 
evolution  for  HO  production  is  also  shown  in  Fig  6.  The  difference 
between  the  SHARC  and  AARC  predictions  can  be  attributed  io  the  use  of  a 
steady  state  solution  to  the  rate  equations  in  AARC.  Ihe  steady  state 
solution  for  HO  requires  the  input  of  s  NO  number  density  profile 
appropriate  for  auroral  conditions.  However,  the  time  dependent  approach 
used  in  SHARC  allows  the  formation  of  NO  consiatent  with  the  chemical 
kinetics  mechanism  describing  the  aurora.  Ther"fore,  SHARC  snd  AARC  are  in 

1J.  O.  Ashihara  and  If.  Takayanagi,  “Velocity  Distribution  if  Ionospheric 
Low-Energy  Electrons,"  PI anet .  Space  Sc l  ,  22,  1201  ( 1 9 7  • ) . 
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NO*  EFFICIENCIES  (IBC  III) 


Q 


Efficiency  (photonc/ i on -po i r) 


figure  '  C'>npm  uon  of  Steady  State  SHARC  and  AARC  Efficiencies  for 

SO ' 

a-t'.at  ant  1  atrrraint  below  120  km  (where  HO  is  in  steady  state),  and  are 
«nitrkii,i t  different  above  120  km  where  SltARC  indicates  that  time  dependent 
effect  a  ate  important  An  illustrative  calculation  of  aurorally  enhanced 
NT  an  1  so*  emission  for  a  vertical  path  to  space  'ros  90  km  is  shown  in 
r  i  f,  *1  The  calculation  weed  an  augmented  HITRAM  data  file,  ^  which  has 
r  I  !  i  *  i ;  ■>  r*.  r  1  lines  for  SO*  and  v  i  hr  a  1 1  ona  1  1  y  excited  HO;  these  lines  arc 
-  i<«-n  from  r  he  data  base  furnished  with  AARC 


t  Rothman,  R.  Samache,  A.  Goldman,  L.  Brown,  R.  Toth,  H. 
Pickett,  R  Poynter,  J.  Flaud,  f.  Camy-Peyret,  A.  Barbe, 

M  Huston,  C.  Fins’ and,  and  M.  Smith,  "The  HITRAN  Database: 
1996  Edition,"  Apn 1  Opt i c  a .  26  6058  (  1989). 
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Figure  8.  Spectral  Radiance  for  a  Class  III  Aurora.  The  Atmospheric 
Path  is  a  Vertical  One  to  Space  from  90  km. 


4 .  CONCLUSION 

SHARC  is  a  new  computer  code  for  calculating  atmospheric  NLTE  radiance 
above  60  km  altitude.  It  is  now  being  distributed  by  GL/OPBj  the  point  of 
contact  is  Dr.  R  Sharma.  Initial  validation  of  SHARC-1  has  been 
accomplished  by  comparison  to  earlier  AFGL  codes,  HAIRM^^  and  AARC.^^  An 
auroral  module,  AURORA,  has  been  incorporated  into  SHARC  and  is  in  the 
process  of  being  tested  and  validated  for  release  at  the  end  of  this  year. 

Additional  upgrades  are  planned.  These  include: 

•  Isotopes  of  the  major  species,  especially  ■‘'■’CO2; 

•  Extension  down  to  50  km  altitude; 

«  Upgrading  the  model  atmospheres; 
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•  Additional  radiators,  e.g.,  OH  and  O2; 

•  A  solar  terminator  model;  and 

•  A  model  for  radiance  fluctuations  due  to  spatial 
variations. 

The  first  three  items  will  be  completed  this  year  and  distributed  with  the 
auroral  option  to  SMARC.  The  last  three  will  be  done  in  subsequent  years. 
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